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ABSTRACT 

We model the mass distribution of long gamma-ray burst (GRB) host galaxies given recent results 
suggesting that GRBs occur in low metallicity environments. By utilizing measurements of the redshift 
evolution of the mass-metallicity (M-Z) relationship for galaxies, along with a sharp host metallicity 
cut-off suggested by Modjaz and collaborators, we estimate an upper limit on the stellar mass of 
a galaxy that can efficiently produce a GRB as a function of redshift. By employing consistent 
abundance indicators, we find that sub-solar metallicity cut-offs effectively limit GRBs to low stellar 
mass spirals and dwarf galaxies at low redshift. At higher redshifts, as the average metallicity of 
galaxies in the Universe falls, the mass range of galaxies capable of hosting a GRB broadens, with an 
upper bound approaching the mass of even the largest spiral galaxies. We compare these predicted 
limits to the growing number of published GRB host masses and find that extremely low metallicity 
cut-offs of 0.1 to 0.5 Zq are effectively ruled out by a large number of intermediate mass galaxies at 
low redshift. A mass function that includes a smooth decrease in the efficiency of producing GRBs 
in galaxies of metallicity above 12+log(O/H)KK04 8.7 can, however, accommodate a majority of the 
measured host galaxy masses. We find that at z ~ 1, the peak in the observed GRB host mass 
distribution is inconsistent with the expected peak in the mass of galaxies harboring most of the 
star formation. This suggests that GRBs are metallicity biased tracers of star formation at low and 
intermediate redshifts, although our model predicts that this bias should disappear at higher redshifts 
due to the evolving metallicity content of the universe. 
Subject headings: Gamma-rays: Bursts: GRB host galaxies 



> 

m 
in 

On 

O 

On 
O 



1. INTRODUCTION 



Th e success of NASA's Swift spacecraft (|Gehrels et al.l 
2004) has resulted in a dramatic increase in the number 
of X-ray localizations of afterglows associated with long 
gamma-ray bursts (GRBs). This increase has resulted in 
a wealth of new information regarding the demograph- 
ics of GRBs and their host galaxies. Investigating the 
environments in which these events occur has long been 
an important path to understanding the nature of GRB 
progenitors, as different origin models have traditionally 
predicted distinct GRB host galaxy populations. The 
connection between GRBs and the death of massive stars 
is now well-established at low redshift (z < 0.3) by the 
association o f GRBs with broad lined SN lc events (for 
a review, see IWooslev fc BloonjlgOOgh. 

Recent observa tions (Ca stro Ceron et all 120081 : 
iSavaelio e t al. 2009|) of X -ray localizations by Swift have 
bolstered previous results showing that GRB host galax- 
ies tend to be bluer, fainter, and more irregular than 
M+ galaxies at similar redshifts dFruchter et alj 119991 
Chary et all l2002t iBloom et all 120021: iLe Floc'h et al I 
20031 iTanvir et all l2004t iFruchter et all 12006 



Castro Ceron et ah! |2006[ ). They tend to have higher 



specific star formation than typical star-forming galaxies 
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(ICharv et alj|2002t iBerger et alJl2003:lChristensen et all 
l2004f) and the peak in their redshift distribution tends 
to broadly track the peak in the overall cosmic star for- 
mation rate of the universe (lBlooml!2003t iFirmani et al.l 



2004 : iNataraian et al 



Kocevski fc Liangl 2006; 



2005 



Jakobsson et al.l 2006; 
Piranl 12001 7 Only 
been associated with grand 
duration GRB has been 



Guetta 



a handful of events have 
design spirals and no long 
associated with an early type galaxy. 

A growing body of spectroscopic evidence has 
also s hown that thes e gala xi es tend to be meta l 
poor dProchaska et all [200l [ Sollcr man et all 1200 
Fruchter et all l2006t iModjaz et all 120061: IStanek et al 
20061: iThone et all 120071 : IWiersema etail 120071 : 
Margutti et al.ll20'07T ) ■ Absorption line spectroscopy has 



in 
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revealed that the regions in which GRB afterglows are 
observed tend to have the low metallicit ies that are 
expected from young stel lar populations (|Fvnbo et al.l 
120031 : Savaglio 'et al.l l2003l). However, there are a few 



exceptions jlFvnbo et al.l 120061: iProchaska et al.l 120071 : 
iFvnbo et all 120081 : IChen et al.ll2008h . The high specific 
star formation rates along with these low metallicities 
are similar to what are seen in star-bursting Lyman 
break galaxies at high redshift. 

There is ample theoretical justification for the a priori 
association of GRBs with short-lived, metal poor pro- 
genitors. The combination of high angular m omentum 
and high stellar mass at the time of collapse (jWoosleyl 
U993tlMacFadven &; Wooslevlll999ft is crucial for produc- 
ing the collimated emission that is required to account 
for the enormous isotropic-equivelent energy released by 
these events. Low metallicity progenitors would, in the- 
ory, retain more of their mass due to smaller line-driven 



2 



D. Kocevski et al. 



stellar winds (|Kudritzki fe Puis! 120001 : IVink fc de Koterl 
I2005T). and hen ce preserve their angular m omentum 
(lYoon fc Langerll2005l: iWooslev & Hegenl200l . 

Recently. iModiaz et alj (2008) showed that a sharp de- 
lineation may exist between the metallicity at the sites of 
broad-lined SN Ic that have been associated with GRBs 
and SN Ic with no detected gamma-ray emission. Us- 
ing a sample of 12 nearby (z < 0.14) broad-lined SN Ic 
without associated GRBs, they found that the chemical 
abundance at the sites of known SN-GRBs (at z < 0.25) 
were systematically lower than those harboring SN with- 
out GRBs, with a boundary between the two samples at 
an oxy gen abundance of roughly 12+log(O/H)KD02 ~ 8.5 
in the iKewlev fc Dopital (|2002f ) scale (see Modjaz et al. 
2008 - Figure 5). This trend is independent of choice of 
the metallicity diagnostic they adopt (see their Figure 6) 
and the mode of SN survey that found the SN without 
GRBs. 

At the same time, the observed trend that many GRB 
host galaxies are less luminous, metal poor, irregular 
dwarf galaxies is in qualitative agreement with the ob- 
served trend of decreasing metallicity of galaxies as a 
function of their stellar mass: the mass-metallicity (M- 
Z) relationship. Alth ough well established at low redshift 
(|Tremonti et aLlfeOO^ ), the M-Z relationship has only re- 
cently been measured for high redshift galaxies where it 
has become clear that the overall normalization of the 
relations hip has decreased throughout the hist ory of the 
universe (|Savaglio et alJ 120051 : lErb et al.ll2"007l . 

As a consequence of the M-Z relationship, any bias 
in the metallicity of the environment that is capable of 
producing a GRB would likely place severe restrictions 
on the type of galaxies that can host such events. While 
earlier studies suggest that the GRB redshift distribution 
tends to broadly track the overall cosmic star formation 
rate of the universe, the question remains as to what 
extent GRB hosts are unbiased tracers of SF in the high 
redshift universe. In this paper, we use empirical models 
based on the measurements of the redshift evolution of 
the M-Z relationship to estimate the upper limit to the 
stellar mass of a galaxy that can harbor a GRB, and 
test the suggestion that GRBs preferentially form in low 
metallicity environments. We detail the prescriptions for 
our model in §2 and expand upon our results in §3. We 
compare our model predictions to published host mass 
values in §4 and discuss the implications of our results in 
§5- 

2. MODEL PRESCRIPTIONS 

To investigate how a potential metallicity cut-off ef- 
fects the resulting GRB host mass distribution we must 
first assume an empirical prescription for the relation- 
ship between a galaxy's stellar mass and its level of 
chemical enrichment. Such a correlation was first ob- 
served by Lequeux et al. (1979); a trend between the 
heavy-element abundance in H II regions and the stel- 
lar mass of irregular and blue compact galaxies. More 
re cently, this corre l ation has been statistically quantified 
bv lTremonti et all (|2004h using of ~ 53,00 galaxies from 
the S loan Digital Sky Survey (SDSS). iTremonti et~aTl 
(|2004h find a tight correlation between galactic stellar 
mass and gas-phase metallicity that spans 3 orders of 
magnitude in stellar mass and a factor of 10 in metal- 
licity. They conclude that the galactic metallicity abun- 



dance rises steeply for stellar masses between 10 8 5 and 
10 10 - 5 M , then flattens for galaxies above 10 10 5 M . 

A basic form of this correlation is a natural conse- 
quence of the conversion of gas to stars within star form- 
ing gala xies, given a mass dependent star form ation ef- 
ficiency (|Schmidtl[l96l ISearle fc Sargentj[l972l ). In the 
context of these simple "closed-box" models, this dispar- 
ity in the efficiency between high and low mass galax- 
ies is thought to be due to variation s in galactic sur- 
face densities as a function of mass dKennicuttl 119981 : 
iMartin fc Kennicuttll2001l : iDalcanton et al.ll2004h . 

It has now become apparent that the effects of su- 
pernovae feedba ck and the infall of metal-poor gas 
( Dalcanton 2007) must also play important roles in shap- 
ing the observed mass-metallicity relationship. Galactic 
winds produced by SNe work to strip galaxies of metal 
enriched gas, with low mass galaxies being more suscep- 
tible to such effects due to their shallower potential wells. 
Energy injection from SNe also heats interstellar gas, de- 
laying the collapse of otherwise cold gas to produce stars. 
At the same time, the infall of metal-poor gas acts to 
dilute the metal content of the ISM. This effect is signifi- 
cant in small galaxies where the infall rate can exceed the 
total star formation rate, causing an net decrease in the 
metallically of the ISM with time. The combined result 
of these mechanisms is that high mass galaxies process 
their primordial gas faster and more efficiently than low 
mass galaxies and are more effective at retaining the re- 
sulting material against wind induced mass loss, leading 
to a positive correlation between stellar mass and metal- 
licity. 

These explanations for the origin of the M-Z rela- 
tionship suggest significant evolution of the relationship 
with redshift. First, one would expect the normaliza- 
tion of the relationship to fall as a function of look- 
back time as metal abundance becomes less common in 
all galaxies. Second, the variations in the efficiency of 
star formation as a function of mass should also change 
the slope of the M-Z relationship as a function of red- 
shift. Efforts to quantify this evolut ion have been the fo- 
cus of sev e ral re c ent observation a l (ISavaglio et al.l 120051 : 
lErb et all l2006t iMaiolino et alj l2008f) and numerical 
(Ide Rossi et alJl2007UKobavashi et al.ll2007l : iTassis et alJ 
2008t iBrooks et al.1 l2007f) investigation. In particular, 
Savaglio et alj ( 2005[ ) used the Gemini Deep Deep Sur- 



vey (GDDS) to examine the M-Z relationship at 0.4 < 
z < 1.0 and found clear evidence for an overall decrease 
in the normalization of the relationship with respect to 
that f ound in the local Universe. Likewise, lErb et aTl 
(2006) utilized 87 rest-frame UV selected star forming 
galaxies to study the nature of the correlation beyond a 
redshift of 2 and came to similar conclusions. 

For our an alysis, we have a dopte d the empirical model 
put forth by ISavaglio et al.l (|2005h to describe the evo- 
lution of the M-Z relationship as a function of redshift. 
This model was developed using their 0.4 < z < 1.0 
G DDS sample al o ng wi th the z ~ 2 galaxies presented 
by IShaplev et a l. ( 2004) to extrapolate the shape of the 
M-Z relationship to highe r redshifts. This emp irical rela- 
tionship (Equation 11 in Savag lio et al.ll2005h allows for 
the average metallicity of a galaxy to be estimated as a 
function of stellar mass at a given redshift and can be 
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Fig. 1. — The e volution of the galaxy mass-metallicity relation- 
ship described by (Savaglio et al. 2005), extrapolated to rcdshifts 
between < z < 5. The overall normalization of the relationship 
is expected to fall with redshift as metal abundance become less 
common in all galaxies. Differential enrichment between low and 
high mass galaxies also leads to an evolution of the relationship's 
slope. The red dotted line represents a low metallicity cutoff of 
12+log(O/H)KK04 = 8-5- Note that our use of Equation 1 beyond 
2 2 is an extrapolation that is beyond the range of the original data 
used to define the model. 

stated as: 

12 + log(O/H) KK0 4 = -7.59 + 2.53 log M* 

-0.097 log M* + 5.17 log t H 
-0.39 log t H - 0.40 log t H log M*, (1) 



where in is the Hubble time and M* is the galactic stellar 
mass. The Savaglio et al. model reproduces several of the 
predicted M-Z relationship properties at high redshift, in- 
cluding the overall reduction in the M-Z relationship nor- 
malization as well as a steeper evolution in the metallic- 
ity of low mass galaxies in comparison to their high mass 
counterparts. Figure 1 shows the metallicity as a func- 
tion of stellar mass for a variety of redshifts as approxi- 
mated by the Savaglio et al. m odel out to z = 5. We note 
that the original data used by 1 Savaglio et all (|2005l ) was 
limited to the range of 8.2 < 12+log(O/H) K K04 < 9.1 
and 0.4 < z < 2.0 and hence the curves at lower metal- 
licities and higher redshifts are extrapolations beyond the 
range of the data used to define the model. 

It is importa nt to examine t he de tails of the diag- 
nostics used by I Savaglio et alj (|2005f ). as different ini- 
tial mass functions (IMFs), for example, can yield fac- 
tor of 2 differences in stellar mass and different metal- 
licity calibrators can likewise result in large discrepan- 
cie s in abundance estim ates. The stellar masses used 
by ISavaglio et al.l (|2005f) to produce their empirical re- 
lationship were estimated through SED modeling of 
multi-band photometry fo r each galaxy, with an initia l 
mass function derived by iBaldrv fe; Glazebrookl (|2003f) . 
Their metallicity values were obtained through nebu- 
lar oxygen abundance estimates calibrated via stellar 
populatio n synthesis and photo i oniza tion models devel- 
oped bv iKobulnickv fc Kewlevl (|2004 hereafter KK04). 
This metallicity diagnostic, which uses traditional strong 
emission li ne ratios, and o t her commonly used cal ibra- 
tions (e.g.. lMcGaughlll991l IKewlev fc Dopital [2001 ) are 
discussed in detail in lKewlev fc Elliso n (2008) who quan- 
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Fig. 2. — The star formation rat e as a function of stella r mass 
between z = — 5 as described by Drory k Alvarez (2008). The 
stellar mass at which this rate turns over evolves smoothly to higher 
masses with increasing redshift. 

tify the systematic offsets amongst the different calibra- 
tions and provide conversion tables. 

We also note that the metallicity value for the bound- 
ary between hosts that harbor broad-lined SN Ic with 
associated GRBs and S N Ic with no d etecte d gamma- 
ray activity reported bv iModjaz et al. (l2008f) was mea- 



sured using the diagnostic proposed bv IKewlev &: Dopital 
(2002) (KD02). In order to co nvert from the KD0 2 
scale to the KK04 scale used by ISavaglio et al.l (|2005h . 
we consulted IKewlev &: Ellison] (|2008[ ) for the appropri- 
ate metallicity calibration conversion (their table 3). We 
find that a value of 12+log(O/H)KD02 = 8.5 approxi- 
mately conve rts to 12+log(G7H)j< K04 ~ 8.66, which we 
quote as the iModiaz et al.l (|2008f ) cutoff metallicity for 
the remainder of the paper. 

In addition to understanding how the average metal- 
licity of a galaxy varies as a function of stellar mass, 
we would also like to know how the number density of 
galaxies and the number of stars being produced in those 
galaxies varies with galactic stellar mass. This will allow 
us to model the effects of a metallicity bias on the overall 
mass distribution of GRB host galaxies, and eventually 
compare those models to the unbiased mass distribution 
of all star- forming galaxies at a given redshift. As with 
the mass-metallicity relation, both the galactic stellar 
mass function and the star formation rate as a function 
of stellar mass are expected to evolve with redshift and 
quantifying this evolution is crucial to understanding the 
distribution of galaxies that are capable of harboring a 
GRB. 

The star formation rate as a function of stellar mass 
(SFRM) in the local universe is well c onstrained. Using 
a sam ple of more than 10 5 galaxies, [Kauffrnau n et all 
(2004) showed that the star formation rate in low mass 
galaxies scales as a power law to their halo mass, peaking 
at roughly log M+ ~ 10.4 Mg, before falling for higher 
mass galaxies. This transition represents the stellar mass 
at which the galaxy distribution changes from younger 
stellar populations and active star forming galaxies to 
systems with older stellar populations and low star for- 
mation activity. 

iDrory k. Alvarez! (120081) used the FORS Deep Field 
survey ( Feulner et al.ll2005h to quantify this relationship 
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Fig. 3. — The galactic mass function as a function of stellar mass 
between 2 = — 5. The number density of galaxies decrease as a 
power law in stellar mass before falling sharply at a characteristic 
mass. The overall number density of high mass galaxies drops 
significantly with redshift. , „ , , 

and its evolution with redshift tor stellar masses and red- 
shifts spanning 9 < log M* < 12 and < z < 5. They 
find that the stellar mass at which the star formation 
rate turns over for high mass galaxies evolves smoothly 
to higher masses with increasing redshift, until the break 
mass disappears entirely and the star formation rate as 
a function of stellar mass can be represented as a sin- 
gle power law. Surprisingly, iDrorv fc Alvarezl (|2008l ) find 
that the power law index representing the low mass re- 
gion of this relationship remains constant even to the 
highest redshifts in their sample. 

For the purposes of this pap er, we have adopted the 
analytic expression presented bv lDrorv fc Alvarea (|2008l ) 
for the star formation rate as a function of stellar mass 
given as: 



M} 



(2) 



where M* represents the break mass at which the star 
formation rate deviates from a power law. We also use 
the best fit parameterizations from Drory et al. (2008) 
for the evolution of the overall normalization and break 
mass with redshift, given as: 

A/o w 3.01(1 + z) 3 - 03 (3) 

Ml «2.7x 10 w (l + z) %1 (4) 

Following IDrorv fc Alvarezl (|2008h . we have fixed the 
power law index to /3 = 0.6 and assume it remains 
constant at all redshifts under consideration. The star 
formation rate as a function of stellar mass between 
< z < 5, as described by Equations (2)- (4), are shown 
in Figure [2] 

The galactic stellar mass function (GSMF) in the lo- 
cal universe is likewise well understood. It has long 
been known that dwarf galaxies represent the largest 
fraction of galaxies in the local universe, with their 
relative number decreasing as a power law with in- 
creasing stellar mass up to some characteristic mass, 
above which the numbe r of galaxies dro ps sharply. At 
low re dshift, the 2dF (ICole et all l200ll) and 2MASS- 
SDSS (|Bell et al.|[200llBlanton et al.ll2003h surveys con- 
strained the parameters of the Schechter function that 
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Fig. 4. — The total star formation rate as a function of stellar 
mass between > z > 5. The portion of the curves highlighted in 
red represent the stellar mass range below the mass limit imposed 
by a metallicity cutoff of 12+log(O/H)KK04 = 8.5. 

is commonly used to describe the distribution of stel- 
lar mass in the Univ erse. The GSMF at high redshift 
has been explored bv[Fbntana et all (l2004h, IDrorv et al.l 
(|2005ft . lConselice et all (|2005ft . and lFontana et all (|2006D 
using a variety of deep surveys, all showing evidence for 
a distinct evolution of the GSMF with cosmic time. Us- 
ing the GOODS- MUSIC catalog of o ver 3000 infrared 
selected galaxies, iFontana et al.l (|2006l ) showed that the 
number density of high mass galaxies drops with redshift, 
while the density of low mass galaxies evolves faster than 
their high mass counterparts out to a redshift of z ~ 1.5. 
The net result of this differential evolution is an increas- 
ing fraction of low mass dwarf galaxies with respect to 
higher mass galaxies at higher redshifts. 

For the purposes of this paper, we have adopted the 
analytic model presented bv lDrorv &: Alvarezl (|2008f ) for 
the galactic stellar mass function given as: 



Af Ti r \ J ]i r ,* ( M \° ( A/\dM 



(5) 



We have again used their best fit parameterizations for 
the evolution of normalization of the mass function as 
well as the characteristic break mass, 



lO: 



4>*{z) w 0.0031(1 + z) 
;M*{z) w 11.35 - 0.22 ln(l 



given as: 

-1.07 



z 



(6) 
(7) 



We further assumed that the power law index below the 
break mass remains constant at a = —1.3 for all redshifts 
under consideration. The GSMF between < z < 5, as 
described by Equations (5)-(7), are shown in Figure [31 

Ultimately, it is important to know the total number of 
stars being produced as a function of stellar mass. Thus, 
we also computed the product of the galactic stellar mass 
function (GSMF) and star formation rate as a function 
of stellar mass (SFRM). This galaxy weighted star for- 
mation rate (WSFR) is shown in Figure 0] at a variety of 
redshifts. The red lines in Figure 2] represent the metal- 
licity biased WSFR, the details of which we will discuss in 
the next section. Between roughly < z < 3, the number 
density of low mass galaxies outweighs that of their more 
massive counterparts, but the cosmic star formation rate 
is largely dominated by these relatively less numerous 
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massive galaxies. The net result is a weighted star for- 
mation rate that peaks at intermediate masses, roughly 
between 10 10 — 1O 11 M . At higher redshifts, the drop 
in the number of high mass galaxies becomes significant 
and the stellar mass function becomes dominated by low 
mass galaxies. At the same time the peak in the SFRM 
decreases smoothly to lower masses with increasing red- 
shift, resulting in a sharp fall in the mass at which the 
weighted star formation rate peaks for z > 3. The mass 
at which the WSFR peaks is plotted as the long dashed 
black line in Figure O If GRBs are unbiased tracers of 
star formation in the universe, and if the y follow the M-Z 
relationship (but see lBrown et al.ll2008f ). then the peak 
of their host mass distribution should roughly follow this 
line. We test this prediction in the following section. 

3. RESULTS 

Using the empirical M-Z relationship expressed in 
Equation 1 , we estimated the stellar mass of a galaxy of 
a given metallicity as a function of look back time. The 
average stellar mass for galaxies with a low oxygen abun- 
dance of roughly 1/3 Zq 1 , or 12+log(O/H)KK04 = 8.5, 
is traced as the red line in Figure O with the red shaded 
region surrounding this line representing the uncertainty 
due to the intrinsic scatter of the M-Z relationship at 
low redshift. Here we have used the values presented by 
Tremonti et al. (2004) to estimate the la scatter about 
the M-Z relationship, and hence the resulting stellar mass 
range, at a redshift of z ~ 0.1. Unfortunately, such de- 
tailed estimates of the scatter associated with the M-Z 
relationship at high redshift are currently lacking and 
therefore for our analysis we assume that this scatter 
is indicative of the scatter at all redshifts under consid- 
eration, which is certainly an over-simplification. The 
region of stellar mass shaded blue and green represent 
the typical masses for dwarf and spiral galaxies respec- 
tively. As expected, the average mass of a galaxy at a 
given metallicity rises as a function of redshift, a direct 
effect of the decreasing normalization of the M-Z rela- 
tionship as a function of look back time. The effects this 
would have on a metallicity biased GRB host popula- 
tion are immediately apparent. If GRBs are limited to 
low metallicity environments, then at low redshift they 
would be relegated to dwarf and low mass spiral galax- 
ies (barring the effects of metallicity gradients, which we 
will discuss in detail in §5"), whereas at high redshift the 
effective mass limit is raised, allowing GRBs to occur in 
a much broa der range of galaxies. A s imilar prediction 
was made bv lRamirez-Ruiz et al.l ()2002l ) based purely on 
theoretical grounds. 

Furthermore, if GRBs are unbiased tracers of star for- 
mation throughout the universe, then their observed host 
mass distribution should cluster about the peak in the 
WSFR represented by the long dashed line in Figure 03 
On the other hand, if they are metallicity biased trac- 
ers of star formation, then their host mass distribution 
should deviate significantly from this curve at low red- 
shifts, peaking instead at the upper mass limit shown in 
red. At some redshift (roughly z ~ 3 for a metallicity 
cutoff of 12+log(O/H)KK04 = 8.5) the stellar mass at 

1 W e assu me a solar metall i city of 12+log(0/H) = 8.7 (|Asplundl 
1200511 in the IPettini &: Page! d2004|) scale and convert tha t value 
to the KK04 scale using Table Thi TKewlev fc Ellison! pOOSft to get 
Zq = 9.0. 
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Fig. 5. — The upper limit on the stellar mass of a GRB host 
galaxy given a sharp metallicity cut-off of 12+log(O/H)KK04 ~ 
8.5 with respect to the typical stellar mass ranges for spiral and 
dwarf galaxies in the local universe. The light red region represents 
the scatter in this limit imposed by the la scatter of the M-Z 
relationship at low redshift. The dashed line represents the stellar 
mass at which the total star formation in the universe peaks at a 
given redshift. Alternatively, the dash-dotted line represents the 
median stellar mass of this distribution, truncated by the upper 
limit set by a metallicity bias. The two regions representing the 
spiral and dwarf galaxy masses overlap between 10 9 — 10 11 Mq. 

which the WSFR peaks crosses this upper mass limit, 
after which the peak in the unbiased and biased mass 
distributions become indistinguishable. 

We also quantified the median of the GRB host mass 
distribution by considering only the product of the 
GSMF and SFRM below the metallicity imposed upper 
mass limit to the GRB host population. The red solid 
lines shown in Figure[4]represent the galaxy masses which 
fall below the upper limit imposed by a metallicity cut- 
off of 12+log(O/H)KK02 = 8.5 for various redshifts. If 
GRBs are metallicity biased tracers of the star forma- 
tion in the universe, then the centroid of this truncated 
total star formation rate would yield the expected me- 
dian stellar mass of a GRB host population as a function 
of redshift. We plot this expected median mass as a dash- 
dotted black line in Figure [5] Although the upper limit 
imposed on the mass of a galaxy that can host a GRB 
increases with redshift, the effects of a galaxy population 
dominated by low mass galaxies along with the shift in 
the type of galaxies producing most of the stars in the 
early universe have the net effect of keeping the median 
GRB host galaxy mass relatively constant with redshift, 
at roughly a mass of 10 8 Mq. Note that this estimate as- 
sumes that the fraction of the total star formation that 
goes into the production of GRB progenitors does not 
change significantly with redshift, host type, or stellar 
mass. This assumption breaks down if environmental 
variables other than metallicity play an important role 
in the formation of a GRB progenitor, normalization 
in the GRB host mass distribution, only their relative 
distribution in stellar galactic mass at a given redshift. 

Unfortunately, the predicted median host mass shown 
by the dash-dotted line if Figure [5] is not currently ob- 
servable, except for low redshift GRBs, which are rare. 
Detection effects and Malquist type biases will lead any 
observational measure of the GRB host mass distribution 
to be biased towards high mass, high surface brightness, 
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galaxies. This would effectively shift the dash-dotted line 
in Figure [5] to higher masses with increasing redshift and 
such completeness considerations have not been incorpo- 
rated into our model. 

4. COMPARISON TO GRB HOST GALAXY OBSERVATIONS 

How do the upper mass limits as inferred from the 
M-Z relationship compare to measured values for the 
subset of the GRBs with known host associations? To 
exami ne this question, we turne d to two recent stud- 
ies by ICastro Ceron. et all (|2008t hereafter CC08) and 
iSavaglio et all (|2009t hereafter SGB09), which compiled 
the galactic stellar masses, star formation rates, and 
dust extinctions for a large sample of GRB host galax- 
ies between < z < 2. CC08 utilized the rest frame 
if -band flux densit ies as i nterpolated from Spitzer's 
(jWerner et all 120041 ) IRAC (|Fazio et all I2004J) and ad- 
ditional NIR observations to obtain an estimate of M* 
for a sample of 30 long-duration GRBs. SGB09 ob- 
tained similar estimates through a combination of optical 
and NIR observations for a sample of 46 host galaxies. 
Both groups used photometric observations in conjunc- 
tion with mass-to-light ratios derived from SED fits to 
measure the total stellar mass of the hosts in their sam- 
ple. The two studies assumed slightly different initial 
mass functions (IMFs) and average mass-to-light ratios, 
introducing a systematic offset between the estimated 
mass values derived from the two samples which we dis- 
cuss in more detail in the next section. 

CC08 and SGB09 found that GRB host galaxies ex- 
hibit a wide range of stellar mass and star forma- 
tion rates, although as a whole they tend towards low 
A/*, relatively dim, high specific star-forming systems , 
confirming previous observations dFruchter et al.l 119991: 



Le Floc'h etaH 120031 : ICharv et alj l2002t iBerger et all 



2003J). TheM* values from the CC08 and SGB09 pa- 



pers are shown in Figures [5] and [7] respectively, with 
upper limits represented by triangular symbols. As in 
Figure 03 the red shaded region in both plots represents 
the upper limit on the stellar mass of a galaxy capable 
of hosting a GRB as imposed by the M-Z relationship 
and its associated la scatter with a metallicity cut-off 
of 12+log(O/H)KK04 = 8.5 It is clear from Figures [6] 
and[7]that a significant fraction of observed host galaxies 
have M+ values that are greater than the predicted upper 
limit to the GRB host mass distribution for such a low 
metallicity cutoff value. Most of these host galaxies can 
be accommodated if the metallicity cut-off is increased 
to 12+log(O/H) KK 04 = 8.7 for the SGB09 sample and 
12+log(0/H) X xo4 = 8.8 for the CC08 sample. Note 
that the resulting spread in the predicted mass range is 
significantly wider for 12+log(O/H)KK04 = 8.8 due to 
the shallower slope of the shallower M-Z relationship at 
this metallicity. Even at this relatively high metallic- 
ity cut-off, with its larger intrinsic spread, three hosts 
in the CC08 sample are still above the predicted mass 
limit, although metallicity gradients within these high- 
mass hosts may explain the existence of these hosts in 
the GRB sample. 

SGB09 and CC08 found median masses of M* ~ 
10 9 3 M Q and M* ~ 10 9 ' 7 M Q respectively, far greater 
than the median host mass predicted by simply look- 
ing at the truncated distribution of total star formation 
as a function of M* (The dash-dotted line in Figure E|) . 
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1 2 3 4 5 

Redshift 

Fig. 6. — Upper limits on the stellar mass of a GRB host galaxy 
given a metallicity cut-off of 12+log(O/H)KK04 = 8-5 ( re d line) 
and 8.7 (blue line) compared to t he masses of 46 host galaxies 
estimated by Savaglio et al. (2009). The dashed line represents 
the stellar mass at which the total star formation in the universe 
peaks at a given redshift. 
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1 2 3 4 5 

Redshift 

Fig. 7. — Upper limits on the stellar mass of a GRB host galaxy 
given a metallicity cut-off of 12+log(O/H)KK04 = 8-5 (ted line) 
and 8.8 (bl ue line) compared to the m asses of 30 host galaxies es- 
timated by C astro Cer o n et al . (2008). The dashed line represents 
the stellar mass at which the total star formation in the universe 
peaks at a given redshift. 

This direct comparison between the expectation peak of 
the WSFR and the median values for the two samples 
is problematic, as detection effects biasing against low 
mass galaxies will heavily influence the observed median 
mass. We can however compare the observed host mass 
distribution to the high end of the expected mass distri- 
bution of all star formating galaxies at a given redshift, 
as detection effects should not effect this region of the 
observed distribution. We address this question in Fig- 
ures M and H where we plot the SGB09 and CC08 host 
mass distributions for galaxies between 0.75 < z < 1.25 
along with the expected unbiased WSFR as a function 
stellar mass at a z = 1 (dashed line). The normaliza- 
tion of the distributions in these plots is arbitrary, with 
the peak of the predicted WSFR and the observed host 
mass distributions both being set to 1. The bracketed 
arrows in Figure [9] represent galaxies for which CC08 
were unable to make firm estimates on the galactic stel- 
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lar mass, resulting only in upper limits accompanied with 
very conservative lower limits. 

It is quite clear that the SGB09 sample is not well de- 
scribed by the expected host mass distribution of unbi- 
ased star forming galaxies at z = 1 . Although the SGB09 
distribution can be expected to artificially fall off at low 
M* due to observational biases, the same cannot be said 
for the lack of high M* galaxies, pointing to an intrinsic 
decline in the GRB host population. The case for the 
CC08 sample is less clear. Their median stellar mass be- 
tween 0.75 < z < 1.25 of = 10 10 - 23 M Q is much more 
consistent with the peak of the unbiased WSFR distri- 
bution, which at z ~ 1 peaks at M* = 1Q 1O - 3O M . This 
median of the CC08 sample does not include the values 
for which only limits exist, which work to broaden the 
distribution to lower M* values, making it less consistent 
with the model distribution. 

We can statistically compare the two observed distri- 
butions to the model distribution by drawing a random 
set of values from the WSFR distribution, equal in size 
to the observed samples, to which we can perform a two- 
sided Kolmogorov-Smirnov analysis. We perform this 
comparison for 1000 trails, using a random realization 
of the WSFR mass distribution in each iteration, while 
measuring the median probability that the model distri- 
bution and the SGB09 and CC08 samples are drawn from 
the same parent populations. For both the SGB09 and 
CC08 samples, the probability that they are randomly 
drawn from the unbiased WSFR distribution is quite low, 
6.3 x 10 -12 and 1.6 x 10~ 05 respectively. Unfortunately, 
the observational biases discussed above leads to the lack 
of completeness at low M+ for both samples making the 
use of a traditional K-S analysis questionable. The me- 
dian WSFR mass will be heavily weighted by low mass 
galaxies, the smallest of which are likely not present in 
the SGB09 and CC08 samples because of detection ef- 
fects. 

At present, without an understanding of the complete- 
ness of the GRB host samples, we can only compare the 
peaks and the high end behavior of the mass distributions 
which we believe should not be effected by observational 
biases. In both cases, the SGB09 and CC08 samples peak 
below the unbaised peak in the galaxy weighted star for- 
mation rate as a function of stellar mass, although the 
discrepancy is much greater when considering the SGB09 
sample. 

5. DISCUSSION 

The comparison between the stellar mass limits im- 
posed by metallicity cut-offs to the measured val- 
ues in the SGB09 and CC08 samples is quite telling. A 
low metallicity cut-off of 12+log(O/H)KK04 = 8.5 is dis- 
favored by current measurements of the stellar masses 
of GRB host galaxies at low and intermediate redshifts. 
However, a comparison of observed GRB host masses still 
appears to favor a metallicity biased mass distribution 
rather than one based solely on the mass distribution of 
all star formation galaxies at similar redshifts. Increasing 
the metallicity cut-off to 12+log(O/H) K K04 ~ 8.7 - 8.8 
allows for the accommodation of most of measured host 
masses, when factoring in the intrinsic spread of the M-Z 
relationship. This is in rough agreem e nt wit h the metal- 
licity cutoff found by iModiaz et~aTl (|2008fl of roughly 
12+log(O/H) KK 04 ~ 8.66 at low redshift (z < 0.25). 



This result is also in gen e ral ag reement with the results 
presented bv iNuza et alJ (|2007h who find a comparable 
metallicity bias through the use of hydrodynamical cos- 
mological simulations in conjunction with assumptions 
of the collapsar event rate. They conclude that the ob- 
served properties of GRB host galaxies are reproduced if 
long GRBs are limited to low-me tallicity progenitors. 

In a similar investigation, IWolf fe Podsiadlowski 
(2007) used the luminosity-metallicity (L—Z) relation for 
galaxies to compare the host galaxy luminosity distribu- 
tions between CC SNe and long GRB host galaxies to the 
expected luminosity function of all star forming galaxies 
at a given redshift. They found that although their ul- 
traviolet based SFR estimates reproduced the CC SNe 
host luminosity distribution extremely well, the same was 
not true for the GRB host population. They found that 
their model SFR estimates would have to exclude lu- 
minous, and hence high metallicity, galaxies in order to 
match the observed GRB host distribution. They con- 
cluded that a metallicity bias with a cutoff of roughly 
12+log(O/H)KK04 ~ 8.7 would be sufficient to repro- 
duce the observed distribution, although they stressed 
that they could not distinguish between a sharp cutoff or 
a decreasing efficiency at producing GRBs as a function 
of decreasing host metallicity with their current data. 
This decreasing efficiency is more realistic than a sharp 
efficiency cutoff and, combined with the spread in the 
M-Z relationship, could explain the existence of outliers 
in Figures [6] and [7l 

Metallicity gradients within galaxies also work to dilute 
observable evidence for a sharp metallicity cut-off in the 
galaxies that can harbor GRBs. The metallicities within 
disk galaxies tend to fall as a f unction of radius from 
the core (e.g. iKewlev et al.1 12005-. and references therein) 
and as such the host integrated light represents an up- 
per limit to the metallicity at the GRB location. The 
nearby galaxies that most closely resemble a typical GRB 
host galaxy at low z for which we have spatially resolved 
spectroscopy are the Large and Small Magellanic Clouds, 
both of which have small inte rnal dispersions of orde r 
0.1 dex in oxygen abundance (|Russell fc Dopital ll992). 
This value is common for star forming dwarf irregular 
galaxies in which metallicity gradients are rather negli- 
gible, although the internal dispersion as measured from 
HII regions in larger g alaxies such as the Milky Way can 
be as high as 0.3 dex (jCarigi et al.l 120051 lEsteban et alJ 
[200l . 

The combined effects of a smooth efficiency cutoff and 
this relatively small expected metallicity gradient on the 
GRB mass distribution are shown as dash-dotted lines 
in Figures [5] and [5J The upper mass limits due to sharp 
metallicity cut-offs are marked by the filled red, green, 
and blue dots as labeled. The dash-dotted lines pro- 
ceeding each limit represents the unbiased WSFR dis- 
tribution convolved with a smoothly broken power law 
decline in the efficiency of producing a GRB at a given 
metallicity, and hence stellar mass. Any effect of a metal- 
licity gradient in a typical host galaxy would work to ex- 
tend the peak of the Z-biased mass distribution to higher 
masses. We assumed that at low M* values, the cutoff 
efficiency is 1, transitioning sharply as M+ — > M cuto ff to 
a power law decline of index a = — 1. We believe that 
such a power-law index can accommodate the spread in 
allowable metallicities from both the effects of a declining 
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Fig. 8. — The GRB host mass distribution as measured by 
SGB09 between 0.75 $J z $C 1.25 compared to the total galaxy 
weighted SFR as a function of galactic stellar mass at z ~ 1 
(dashed line). The mass limits due to sharp metallicity cut-offs of 
12+log(O/H) KK 04 values of 12+log(O/H) KK 04 = 8.5, 8.6, and 8.7 
are represented by a red, green, and blue filled circles respectively. 
The combined effects of a smooth efficiency cutoff in the produc- 
tion of a GRB as a function of metallicity are shown as dash-dotted 
lines proceeding each upper mass limit. The peak of the SGB09 
sample is roughly an order of magnitude below the expected peak 
of an unbiased galaxy weighted SFR. 
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Fig. 9. — The GRB host mass distribution as measured by 
CC08 between 0.75 ^ z ^ 1.25 compared to the total galaxy 
weighted SFR as a function of galactic stellar mass at z ~ 1 
(dashed line). The mass limits due to sharp metallicity cut-offs 
of 12+log(O/H) KK 04 values of 12+log(O/H) K K04 = 8.5, 8.7, and 
8.8 are represented by a red, green, and blue filled circles respec- 
tively. The arrows represent galaxies for which CC08 could only 
estimate conservative upper and lower limits to their mass. The 
combined effects of a smooth efficiency cutoff in the production of 
a GRB as a function of metallicity are shown as dash-dotted lines 
proceeding each upper mass limit. The CC08 mass distribution is 
much broader than the SGB09 sample at this redshift, although the 
galaxies for which only upper and lower limits exists pull the peak 
of their distribution between 10 9 — 10 10 Mq, below the expected 
peak of an unbiased galaxy weighted SFR. 



efficiency and the small metallicity dispersion expected 
in GRB host galaxies. 

In the context of these two effects, the resulting M cuto ff 
now can be understood as the peak in the predicted GRB 
host mass distribution at low redshift and not a sharp 
upper limit. As such, this smooth decrease in efficiency 
can accommodate host galaxies of much higher stellar 



mass than the scatter in the M-Z relationship alone. At 
a metallicity cut-off of 12+log(O/H)KK04 = 8.7, for ex- 
ample, galaxies of M* ~ 10 11 M Q are permitted by the 
model (in relative abundance) whereas the scatter in the 
M-Z relationship with a sharp efficiency cut-off would 
strictly exclude galaxies above M* ~ 10 10 at a z ~ 1. 

A cut-off of 12+log(O/H) KK 04 - 8.7, or 1/2 Z Q , does 
contain most of the CC08 sample, although the low mass 
location of the peak in the SGB09 sample points to a 
systematic difference between the two samples. The two 
host samples have a total of 25 overlapping objects and 
CC08 discussed the differences between the two studies 
in some detail. They concluded that the higher median 
mass in their sample reflects a lower mass-to-light ratio 
obtained from the subset of galaxies for which they per- 
formed SED fits, compared to the average value obtained 
by SGB09 through SED fitting to their entire host sam- 
ple. They also added that the use of optical-NIR SEDs by 
SGB09 may underestimate the effects of dust extinction 
for obscured galaxies, giving rise to further discrepancies. 

We also note tha t CC08 used a traditional Salpeter 
IMF (|Salpeterlll955l ) to estimate the relative number of 
low mass, and hence undetected, stars within a galaxy, 
whereas SGB09 utilized a modified S alpeter IMF pre- 
sented bv lBaldrv fc Glazebrookl (|2003l ) . The host masses 
derived through the use of these two IMFs can differ sub- 
stantially, as the traditional Salpeter IMF tends to over- 
estimate the number of low-mass dwarf stars compared 
to up d ated models p resent ed by Baldrv &: Glazebrookl 
(|200a . lKroupal (|200ll) . and iBochanskil (|2008fl . We esti- 
mated that masses derived through the use of a Baldry & 
Glazebrook IMF are systematically lower by roughly 85% 
compared to those found through the use of a Salpeter 
IMF for a given mass-to-light ratio. This combined with 
the lower effective mass-to-light ratio used by SGB09 
may explain the discrepanc ies between the two s amples. 
It is important to note that I Savaglio et alj (|2005f ) explic- 
itly used the Baldry & Glazebrook IMF to obtain the M- 
Z relationship parameterization that we use in Equation 
[TJ therefore their sample makes for a more meaningful 
comparisons to our model. 

The effects of a metallicity bias in the GRB progenitor 
population combined with an evolving M-Z relationship 
would suggest that afterglow associations with massive 
galaxies of > 10 11 M Q should be limited to high red- 
shifts events. Unfortunately, very few host galaxies have 
measured masses above z > 2 to test this directly, despite 
the median redshift of Swift detected GRBs of z ~ 2.75, 
highlighting the difficulty in observing many of these high 
redshift hosts. Despite this increase with redshift in the 
upper limit in the mass of galaxies capable of harboring 
a GRB, we find that an evolving galaxy populations in 
which dwarf galaxies represent a larger fraction of the 
star forming galaxies in the distance universe results in 
a median GRB host mass that remains fairly constant 
between < z < 3 Above z > 3, we see this upper limit 
fall sharply, implying that a large fraction of GRBs at 
high redshift will still occur in low mass galaxies. 

Finally, if the normalization of the M-Z relationship 
for galaxies decreases as a function of lookback time, 
then there should be some redshift at which a metal- 
licity biased galaxy populations would become indistin- 
guishable from the star-forming field galaxy population. 
For a metallicity cutoff of 12+log(O/H)KK04 = 8.7, we 
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find that the peak in the stellar mass distributions be- 
tween these two populations should become equal at a 
rcdshift of z ~ 2, with the biased and unbiased popu- 
lations becoming less distinguishable at higher redshifts. 
The greatest discrepancy between a metallicity biased 
host population at that of the population of all star form- 
ing field galaxies would occur at low to intermediate red- 
shifts. This may explain the discrep ancy between high 
reds hift host properties reported by iChen et all (2008) 
and iFvnbo et alj (120081) and the properties reported by 
iWolf &; Po dsiadlowski (2001]) for hosts at intermediate 
redshifts. IChen et al.1 <(2008fi found that the UV lumi- 
nosity distribution of long GRB hosts is largely con- 
sistent with their being drawn from a UV luminosity 
weighted random g alaxy population at similar redshifts. 
IFvnbo et al.l (|2008l ) reported on similar agreements when 
comparing the luminosity and metallicity distributions 
of GRB hosts to UV-selected star forming galaxies at 
z ~ 3. This is i n stark disagreement with th e conclu- 
sions reported by I Wolf &: Podsiadlowskil (|2007l ) who find 
that a metallicity truncated field population is required 
to match the luminosity distribution of GRB host galax- 
ies at redshifts of 0.4 < z < 1.0. This dichotomy between 
the high redshift and low redshift comparisons would be 
expected, if at some point, the two populations become 
indistinguishable as the average metallicity of the field 
galaxies falls with increasing lookback time. 

6. CONCLUSIONS 

We find that dearth of massive GRB host galaxies 
at low and intermediate redshifts exceeds that expected 
from the decline in the predicted number of massive star 
forming galaxies at similar redshifts. We therefore con- 
clude that there is sufficient evidence to indicate that 
GRB host galaxies are metallically biased tracers of star 
formation at low and intermediate redshifts and suggest 
that this bias should disappear at higher redshifts due to 



the evolving metallicity content of the early universe. We 
find that a galaxy mass function that includes a smooth 
decrease in the efficiency of producing GRBs in galax- 
ies of metallicity above 12+log(O/H)KK04 = 8.7 accom- 
modates a majority of the measured host masses. This 
is in rough agreement with the metallicity cutoff found 
bv lModiaz et all (|2008f) of roughly 12+log(O/H) KK 04 ~ 
8.66 at low redshift (z < 0.25). Throughout our analy- 
sis, the modeling and subsequent metallicity comparisons 
have been performed in the same, consistent fashion and 
in the same metallicity calibration scale, in order to avoid 
any systematic differences between the various metallic- 
ity diagnostics used in the literature. 

For a metallicity cut-off of 12+log(O/H)KK04 ~ 8.7, 
the predicted peak in the GRB host mass distribution 
and the stellar mass at which the weighted star formation 
rate peaks become equal at z ~ 2, with higher values of 
12+log(0/H)cutofi pushing this intersection to lower red- 
shift. This limits the redshift range in which the differ- 
ences between a metallicity biased GRB host population 
and that of unbiased star forming galaxies can be tested 
through direct luminosity or mass distribution compar- 
isons. Therefore, comparisons of these distributions at 
low and intermediate redshifts will be crucial to further 
inquires into the nature of the metallicity bias in the 
GRB host population. 
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